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We discuss the connection between the contributions of large field configurations and the large 
order behavior of perturbation theory. For quenched QCD, the sensitivity of the average plaquette 
to a removal of large field configurations has a narrow peak near j3 =5.6. Various analysis of 
the order 10 weak coupling series for the plaquette give robust indications for a singularity in the 
third derivative of the free energy (second derivative of the plaquette) with respect to j3, near j3 
= 5.7. We report results of numerical calculations. The peak in the third derivative of the free 
energy present on 4 4 lattices disappears if the size of the lattice is increased isotropically up to a 
10 4 lattice. On the other hand, on 4 x L 3 lattices, it persists when L increases. The location of the 
peak coincides with the onset of a non-zero average for the Polyakov loop and seems related to 
the finite temperature transition. We also discuss the discrepancy between the perturbative series 
and the numerical values of the plaquette. 
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1. Introduction 

Perturbation theory can be a frustrating tool for field theorists. Sometimes, it provides ex- 
tremely accurate answers, sometimes it is not even qualitatively correct. In recent years, our main 
goal has been to construct modified perturbative series which are converging and accurate. As 
briefly reviewed in Section |2], our approach consists in removing large field configurations in a way 
that preserves the closeness to the correct answer. 

In the case of quenched QCD, there are several questions that are relevant for this approach 
and that have been addressed. How sensitive is the average plaquette P to a large field cutoff |jl|]? 
How does P behave when the coupling becomes negative [^]? How does P differ from its weak 
coupling expansion [f3|, ^|]? Are all the derivatives of P with respect to /3 continuous in the crossover 
region? The analysis [Q, of the weak series for P up to order 10 [||] suggests an (unexpected) 
singularity in the second derivative of P, or in other words in the third derivative of the free energy. 
In the following, we report our recent attempts to find this singularity. As all the technical details 
regarding this question have just appeared in a preprint [Q] , we will only summarize the main results 
leaving room for more discussion regarding the difference between series and the numerical values 
of P. 



2. Large field configurations and perturbation theory 

The reason why perturbation theory sometimes fail is well understood for scalar field theory. 
Large field configurations have little effect on commonly used observables but are important for the 
average of large powers of the field and dominate the large order behavior of perturbative series. A 
simple way to remove the large field configurations consists in restricting the range of integration 
for the scalar fields. 
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For a generic observable Obs. in a A 4 theory, we have then 

K 
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The method produces series which apparently converge in nontrivial cases such as the anharmonic 
oscillator and D = 3 Dyson hierarchical model ||]. 

The modified theory with a field cutoff differs from the original theory. Fortunately, it seems 
possible, for a fixed order in perturbation theory, to adjust the field cutoff to an optimal value 
tymax (A , K) in order to minimize or eliminate the discrepancy with the (usually unknown) correct 
value of the observable in the original theory. In a simple example [g], the strong coupling can be 
used to calculate approximately this optimal <p limx (X,K). This method provides an approximate 
treatment of the weak to strong coupling crossover and we hope it can be extended to gauge theory 
where this crossover JT0| ] is a difficult problem. The calculation of the modified coefficients remains 
a challenge, however approximately universal features of the transition between the small and 



large field cutoff limits for the modified coefficients of the anharmonic oscillator [11], suggest 
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the existence of simple analytical formulas to describe the field cutoff dependence of large orders 
coefficients. 

This method needs to be extended to the case of lattice gauge theories. Important differences 
with the scalar case need to be understood. For compact groups such as SU (N), the gauge fields are 
not arbitrarily large. Consequently, it is possible to define a sensible theory at negative j3 = 2N/g 2 . 
However, the average plaquette tends to two different values in the two limits g 2 —> ±0 [g]. This 
precludes the existence of a regular perturbative series about g 2 = 0. A first order phase transition 
near j3 = —22, was also observed [Q] for SU (3). 

The impossibility of having a convergent perturbative series about g 2 = is well understood 
fll2[ ] in the case of the partition function for a single plaquette which after gauge fixing to the 
identity on three links reads. 

Z= [ dUe-W-z ReTrU \ (2.1) 



If we expand the group element U = e' gA with A=A a T a and the Haar measure in powers of g, we 
obtain a converging sum that allows us to calculate Z accurately, however, the "coefficients" are g- 
dependent. This comes from the finite bounds of integration of the gauge fields that are proportional 
to l/g. If g is small and positive, we can extend the range of integration to infinity with errors 
that seem controlled by e -2 ' 3 . By "decompactifying" the gauge fields, we have transformed a 
converging sum into a power series in g with constant coefficients growing factorially with the 
order. The situation is now resemblant to the scalar case and can be treated using this analogy. 
We can introduce a gauge invariant field cutoff that is treated as a g independent quantity. For a 
given order in g, one can use the strong coupling expansion to determine the optimal value of this 
cutoff. This provides a significant improvement in regions where neither weak or strong coupling 



is adequate [12]. 

This program can in principle be extended to LGT on D-dimensional lattices, however the cal- 
culation of the modified coefficients is difficult. An appropriately modified version of the stochastic 
method seems to be the most promising for this task. As the technology for completing this task is 
being developed, we will discuss several questions about the average plaquette and its perturbative 
expansion. 

3. The average plaquette and its perturbative expansion in quenched QCD 

We now consider a SU (3) lattice gauge theory in 4 dimensions without quarks (quenched 
QCD). We use the Wilson action without improvement. Our main object will be the average 
plaquette action denoted P and can be expressed as — d(ln(Z)/6L 4 )/dj6. The effect of a gauge 
invariant field cutoff is very small but of a different size below, near or above j3 = 5.6 (see Fig. 6 
of Ref. [|I|]). This is in agreement with the idea that modifying the weight of the large field config- 



urations affects the crossover behavior [13]. The weak coupling series for P has been calculated up 
to order 10 in Ref. [g]: 

10 
m=l 

The coefficients are given in table 1 . The values corresponding to the series and the numerical data 
calculated on a 16 4 lattice is shown in Fig. ||. A discrepancy becomes visible below j3 = 6. The 
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WEAK ORDER 10; PADE 4 over 6 



STRONG PADE 7 over 7 





Figure 1: Regular weak series (blue) and 4/6 weak Pade (red) for the plaquette (left); 7/7 strong Pade (right) 



FIT ORDER 8 




FIT ORDER 10 




6.2 6.4 



Figure 2: Logw\P — Pw\ for order 8 (left) and 10 (right, in a different range of j3); the constant is fitted 
asumming a 2 (blue) or a 4 (red). 



situation can be improved by using Pade approximants, however, they do not show any change in 
curvature and often have poles near j8 = 5.2. For comparison, Pade approximant for the strong 



coupling expansion [ 14 ] depart visibly from the numerical values when /3 becomes slightly larger 
than 5. In conclusion, it is not clear that by combining the two series we can get a complete 
information regarding the crossover behavior. 

The difference between the weak coupling expansion ify and the numerical data P can be 



further analyzed. From the example of the one-plaquette model [12], one could infer that by adding 
the tails of integration, we should make errors of order e~ c ^, for some constant C. Consistently 
with this argument, the difference should scale as a power of the lattice spacing, namely 



PNonPert. = (P ~ Pw) x <^ oc ( e 



33 P 



(3.1) 



A case for A = 2 has been made in Ref. [^] based on a series of order 8. Another analysis supports 
A = 4 (the canonical dimension of F^yF^) [Q, [l5|]. Fig. || shows fits at different orders and in 
different regions that support each of these possibilities. It would be interesting to study cases 
where long series are available and non-perturbative effects well understood in order to define a 
prescription to extract the power properly. 

The series Pw has another intriguing feature: r m = b m /b m -\, the ratio of two successive co- 
efficients seem to extrapolates near 6 when m — > °° when m becomes large [ffl]. This suggests a 
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behavior of the form 

P=(1//3 C -1//3)^(A + A 1 (/3 C -/3) A + ....), 

as encountered in the study of the critical behavior of spin models. We have reanalyzed [§] the 
series using estimators JI^ ] known as the the extrapolated ratio (R m ) and the extrapolated slope 
(S m ) in order to estimate j3 f and 7. We found that the weak series suggests 

Poc (1/5.74- l/j3) 108 . (3.2) 



These estimators are sensitive to small variations in the coefficients and show a remarkable stability 
when the volume is increased from 8 4 to 24 4 . The numbers are in good agreement with the esti- 
mates of Ref. [Q] with other methods. A finite radius of convergence is not expected and one does 
not expect any singularity between the limits where confinement and asymptotic freedom hold. It 
may simply be that the series is too short to draw conclusion about its asymptotic behavior. A 
simple example where this happens [Bl] is 



QiP) 



dte- , t a [l-tp c /(ap)]-r 



(3.3) 



with a sufficiently large. If m << a, r m ~ j8 f (l + (7 — l)/m), For m » a we have r m <=c m and 
the coefficients grow factorially. 

If we take Eq. (3.2) seriously, it implies that the second derivative of P diverges near /3 =5.7. 
We have searched for such a singularity [§J. We have shown that the peak in the third derivative of 
the free energy present on 4 4 lattices disappears if the size of the lattice is increased isotropically 
up to a 10 4 lattice. On the other hand, on 4 x L 3 lattices, a jump in the third derivative persists 
when L increases. Its location coincides with the onset of a non-zero average for the Polyakov 
loop and seems consequently related to the finite temperature transition. It should be noted that 
the possibility of a third-order phase transition has been discussed for effective theories of the 



Polyakov's loop [17]. 



A few words about the tadpole improvement [18] for the weak series. If we consider the 
resummation 

K 

-K-\\ 



P w (\/P)= 1 £e m f3 R m + 0(P 



R 



(3.4) 



m=l 



with j6# = j3(l — Y,m=\b m fi m ), the ratios e m /e m -\ stay close to -1.5 for m up to 7, but seem to 
start oscillating more for large m. 



m 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




2 


1.2208 


2.9621 


9.417 


34.39 


136.8 


577.4 


2545 


11590 


54160 


&m 


2 


-2.779 


3.637 


-3.961 


4.766 


-3.881 


6.822 


-1.771 


17.50 


48.08 



Table 1: b m : regular coefficients; <?,„: tadpole improved coefficients 



This research was supported in part by the Department of Energy under Contract No. 
FG02-91ER40664. We thank G. Burgio, F. di Renzo and P. Rakow for interesting discussions. 



258/5 



Third-order phase transition in QCD? 



Y. Meurice 



References 



[1] L. Li and Y. Meurice, Effects of large field cutoffs in scalar and gauge models, Nucl. Phys. Proc. 



Suppl. 140 (2005) 788-790, Qhep-lat 704 90 9 60 - 
[2] L. Li and Y. Meurice, Lattice gluodynamics at negative g**2, Phys. Rev. D 71 (2005) 016008, 



Dhep-lat/0410029[ l- 

[3] G. Burgio, F. Di Renzo, G. Marchesini, and E. Onofri, Lambda**2-contribution to the condensate in 



lattice gauge theory, Phys. Lett. B422 (1998) 219-226, [hep-ph/9706209]. 



[4] R. Horsley, P. E. L. Rakow, and G. Schierholz, Separating perturbative and non-perturbative 



contributions to the plaquette, Nucl. Phys. Proc. Suppl. 106 (2002) 870-872, [ |hep-lat/0110210 1 
[5] L. Li and Y. Meurice, About a possible 3rd order phase transition at t = in 4d gluodynamics, 



hep-lat/0507034 



[6] F. Di Renzo and L. Scorzato, A consistency check for renormalons in lattice gauge theory: 



beta**(-10) contributions to the su(3) plaquette, JHEP 10 (2001) 038, [hep-lat/0011067]. 



[7] Y. Meurice, A simple method to make asymptotic series of Feynman diagrams converge, Phys. Rev. 



Lett. 88 (2002) 141601, [tiep-th/0103134|. 



[8] L. Li and Y. Meurice, A tractable example of perturbation theory with afield cutoff: The anharmonic 



oscillator, J. Phys. A 38 (2005) 8139-8153, flhep-th/0506038| ]. 
[9] B. Kessler, L. Li, and Y. Meurice, New optimization methods for converging perturbative series with a 



field cutoff, Phys. Rev. D69 (2004) 045014, flhep-th/0309022Q 



[10] J. B. Kogut and J. Shigemitsu, Crossover from weak to strong coupling in su(n) lattice gauge theories, 
Phys. Rev. Lett. 45 (1980) 410. 

[11] L. Li and Y. Meurice, Asymptotically universal crossover in perturbation theory with afield cutoff, 



hep-th/050719( 



[12] L. Li and Y. Meurice, An example of optimal field cut in lattice gauge perturbation theory, Phys. Rev. 



D71 (2005) 054509, [hep-lat/0501023] 



[13] G. Mack and V. B. Petkova, Comparison of lattice gauge theories with gauge groups z(2) and su(2), 
Ann. Phys. 123 (1979) 442; G. Bhanot and M. Creutz, Variant actions and phase structure in lattice 
gauge theory, Phys. Rev. D24 (1981) 3212; P. de Forcrand and O. Jahn, Comparison of so(3) and 



su(2) lattice gauge theory, Nucl. Phys. B651 (2003) 125-142, [hep-lat/0211004]; 



M. Hasenbusch and S. Necco, Su(3) lattice gauge theory with a mixed fundamental and adjoint 



plaquette action: Lattice artefacts, JHEP 08 (2004) 005, [hep-lat/0405012] 



[14] R. Balian, J. M. Drouffe, and C. Itzykson, Gauge fields on a lattice. 3. strong coupling expansions and 
transition points, Phys. Rev. D19 (1979) 2514. 

[15] P. Rakow et al. these proceedings. 

[16] B. Nickel, Phase Transitions, Cargese 1980. Plenum Press, New York, 1982. 

[17] A. Dumitru, J. Lenaghan and R. Pisarski, Deconfinement in matrix models about the Gross-Witten 



point, Phys. Rev. D71 (2005) 074004, [ |hep-ph/0410294[ l. 
[18] G. P. Lepage and P. B. Mackenzie, On the viability of lattice perturbation theory, Phys. Rev. D48 



(1993) 2250-2264, [hep-lat/9209022| 



258/6 



